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o  role   of  the    suberitical  assembly  is  gradually  being 
enlarged  to  include  studies  of  transient  or  lc   responses. 

Demands  for  an   inexpensive  and   qui  thod  for  c. 

the  effect  of  various   parameters  upon  control  response  have 
led  I  rfclonn   into  the  utilization  of   too   a*  iwbly  for 

such  pur*>osos.      Nh  the  most  important  adv  0   of  the 

suberitical  assembly  is   the   eane  with  which  expei'iments  may 
be  performed.  ../sical   r-afeguarda  and  restrictions   are 

at   a  minimum,    and  no  extensive  shielding   Is   re  .        ;ie 

relatively  small  sise  of  the   suboritical  assembly  also  tends 

It  easier  and  less  expensive   to  vary  the  desired 
parameters. 

rrovlng  fron  s srvwiu otoan I—  theory  and  treating  the 
assembly  aa   a  "black  box",   many  investigators  have   aptly 
applied  the  concept   of  the  transfer  function  to  describe  the 
response  to   a  given  forcing  function.      In  practice,   phase   and 
Samplitude  relations,    uhon  cotisidered  as  functions   of  the 
driving  fre-uency,    are  collectively  designated  as   th«   as.- 
tj'anafer   function. 

The  parametric  forcing  or  driving  function  is  of  two 
types,    reactivity  forcing   and  source  forcing.      In  general, 


reactivity  faro  Log  is  a  In  neutron 

lsity   c!ue   to  a  nraa  >nic  variation  5  -blvity, 

.ry  acco-  - 
pllshed  by  oscillating  sons    I  o© 

forcing,   in  conoraat  to   reactivity  forcing,    In  no  wa^  cts 

or  change*  the  asoosibiy  ft  y.      It  is  nil  ed 

I  iatiOii  j. 

The  objective  of  this   thesis  is  to  present  the  sub- 
critioal  assembly  trans for  function  which  was  rrained  by 

I   technique  of  source  f 


tor-y  of   the   suborltioal  pile   and  Its  role   in   the 

davalc  :   of  the   first  self -sustaining  chain  reaction  has 

beoone    legend   (1),   taut  work   Is    cotitlnuin  ike    the 

a»r  vehicle   m  orely  an  object  of  In- 

ves  its©Xf» 

The  initial   servo  jjlsn  concept  was   presented  by 

,  p    and  Owens   (9   and  10)   at  the   Knol  ic 

.  atory.  I   the  roots  of  the  kinetics  equation 

are  negative,   the  nfer  function  was  noted  as   oeing 

•  sane  as  that  for  a  BlAlwon  phase  shift  network.      This 

permitted  It  to  be  represented  by  a  ran.  !• 

■ 
kinetic   similator,   whose  voltage  respon-e   represented  the 

le  neutron  flux  respo ■■.■    ,   was  succec  .y  utilised  for 

■   tert  to  doterraine  the  stability  of  an  automatic  control 

system   and   Its    ability   to  maintain  a  desired  power   level. 

Franz    (5)   determined  substantial!  same   type   t 

far   fttnstloi         ft  used   a  move  sophisticated  approach. 

definition,    the  tvi  relates   sinusoidal  varia- 

tionr  it  can    >e  handled  standard  techniques   of 

■  ■IVOHMOl  L.heory.      The  ratio  of  tr.n   Laplace   trail-  of 

the  output   to  that  of   t  ';  was  found   to  give  the  ;fei 

function  wore   convenient ly.      In  both   the    above  si  Ions   it 


was  oarofu1  it  the  functions  dU     Lned  wore  thoss 

of  a  pile  simulator  and  described  t;       responses  only  to 

I  same  degree  that  the  equations  Utilised  described  the 
conditions  in  the  pile. 

It  remained  for  ilarrer,  .oyer,  aru.  ;.rucoff  (6),  based  on 
their  now  classical  experiment  on  the  Argot     r-2  reactor, 
to  justify  the  uae  of  the  transfer  function  to  describe  a  re- 
actor as  a  control  loon  element,   Since  their  experiment  it 
has  become  start-     rocedure  to  sfer  func- 

tion of  critical  ■  ctors  for  stability  cri- 

teria. This  method  lias  also  been  applied  to  I       of  pro- 
sed reactors,   volar  .,   dth#  and  Kico  (3)  found  the  aero 
ansfer  function  of  the  1    -II3  ■  c  u  assembly  of 
-  ,  a  fast  critical  assembly.  The  results  were  U3ed  to 
determine  the  best  oscillator  rods  and  measurement  techniques 
to  ne  used  on  the  actual  as;     ,  when  constructed. 

In  all  the  references  found,  the  ir     tlon      rily 
pertained  to  a  critical  assembly  and  to  the  solution  of  the 
transfer  function  by  reactivity  fore  in;-..   I  -,ar   than  a  cur- 
sory discussion  in  Glasstono  (6,  p.  2>2) ,   no  unclassified 
reference  was  found  which  contained  any  information  con- 
cerning the  transfer  function  of  a  subcrltical  assembly  as 
such  or  concerning  source  forcing. 
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The  response  of  the  sukoritioal  ub  socio  ly  to  source 
forolng  oan  ->e  derived  from  the  tlta* -do  pendent  general  diffu- 
sion e  uattion.      /or  simplicity  all  neutrons  are  oonslderod 
therwelt   and  hence  the  equation  can  be  written  as  (7»      • 
101) 

B  *,2  0  -    2     0  *   3  *  -T?     a  ~*r  -T*  .1 


^-F*  «t»*»w-3|     -™  c)t 

vbe^'O  at  Is  the  mean  velocit  'she  t  I  neutrons •     In 

general  the  symbols  In  the  derivation  are  In  c         i  usage  (6, 
pp.  .       at  the  i  ,g*  are  given  of 

S^fwbols  for  .lent  refer. 

The  eourcu  term  is  m&d*  ,       '.*e  firat 

portion  is  the  rom  fission. 

to  th*s  i%rai  age  theory  this  prompt  neutron  source 
torn  may  be  expressed  as 

S  Ml*  (3  I      «    2^*  e*"     T  JSq,  2 

where   ft  is  the  tot a  irons  in  any 

one  generation.  Batons'  li   fefe*  relayed  neutrons 

which  resu  fission  ;  the 

deoay  of  certain  fission  products*     Generally*  six  a  are 

reeogr  .  eaoL  .  -  a  disoret  lay  MttaiMttt    A, 
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Lag  cr.."* 

.otion  of  average  |     of  delayed  neutron* 

Fractl  n  of  l6*1       of  delayed  neutrons 

Concentration  .f  average  group  of  do 
neutron  pi -eoursort 

oentratlon  of  !***  :rou;i  of  delayed  neutron 

- 

Thermal  neutron  diffusion  coefficient  c«« 

kff active  Rulttrjlicition  factor 
Multiplication  factor  for   an  Infinite  assembly 
Seutron  diffusion  length  in  the  moderator         cm* 
average  neutron   llfet.         in  a  finite  assembly     sec* 
Mean  neutron  lifeline   in  an  Infinite  assembly     sec. 


Decay  constant  of  av  group  of  delayed 

neu t  ren  |  i*e cur  aor  s 

Decay   com  t  ant  of  the  l**3  group  of  delayed 
neutron  precursors 

lfc»b©r  of  neutrons  of  thermal  energies 

Resonance  escape  probability 

iherraai  neutron  flux 

t«  of  production  of  source  neutrons 

Macroscopic  absorptl  n  crosssootlon 

Titne 

Fermi  ag* 

Thermal  neutron  mean  velocity 

»oncy  of  oscillation 


see 


•1 


seo« 


•m  \ 


neutrons 


om< 


seo. 


neutrons 


0      . 


-1 


a«?o. 
cm.** 

©isu/s 

radia 
sec. 


•ro- 


X 


T 


6 

fract:   (3   (11*   .    ).       C*  defined  as      onc©n- 
.tion  of  the  i^b  ^roup  of  utron  emit:    ,  ;  M 

rat©  of  formation  of  the  precursor  i  on  by 

e  rat©  of      ctior.  of  the  delayed  neutrons  is 

BOP'     lly  equal  to  the  decay  rate  of       o cursor,  this 

© 

rate   equals     ^L.    X,    ,    for  the  six  groups  aro&ei    considera- 

tion.  ermi   tig©   of   a  delayed  neutron  ass^  *al 

to    th.it   of  a  t  neutron,    the   rlowi-    -      .n   nonloakag© 

probability  will  b©  e  .     This  is  not  strictly  tru<-  ..use 

the  l :  gy  of  the  delaycsd  neutrons    f  orally  less  she 

average   energy  of  the   prompt  neutrons*   but   the   approximation 
is  cons id ©rod  reasonable    (6,      .        /).         m   delayed  neutron 

te  term  nay  be   written  as 

© 

3d  *  -f  ^lXi:i  ■*   k 

where  -jo    U         Ined  a  •  reac  I  •»••  L  iity. 

third  portion  is   th«  external  >  term  which  will  oe 

defined  as 

ju>t 

Bq.   5 

where      ,      La  orcp lex  amplitude. 

Per   further  simplification  of  the  natuenatlcs,    the   s^ 
and   tlta©   va_  es  will  be  considered  separable  so   that   the 

partial  differential  equations  -uced  to  ordinary 


>    T 


differential  equations    (6,   f  .  )• 

With  substi      t:  ^   if  ,ionfi   i  ,    ...,    and  $  i::tb  ^uatlon 

',  lulting  egression  is 


first  M  1  bo  considered  for   t  i  d  solu- 

tion,  and  hence   V%m  -  fy(#j   p.   33).  ill  in    jt     U 

mean  lifetime  of  tie    t  .1  neutrons  In  an  Infinite 

medi      ,  it  equals  *~ .  division  by    2.     m 

the  above   substitutions,  -.tion  6  beiiwiwi 


e  e  -t  jp  -     X    ^ 

.{\rir+\\t  e  (l-  (jjk^e         0  +  ~f~-  i  2.  A 


S  °   dt 

*■  a 

let  ,     e'     T 

k  ■  J? „ 

1  ♦ 
and  B  I 

X .  -^rj. ,  n.  v 

1  +    1,1 
where  k  ti   the  effoetive  nuiti  .tor   end  «*    is   the 

av  arm&l  m  Laie,  both  in  a  finite  aaaenoly. 

d  of  k  and   JL  $ 
it   is  possible  to  writ*  ?  aa 


J*    «.-.i» 


ItiedUB   : 


M*  3. 


•                               i<tft 
-  ^  4  (i  -(3)k#  *  -4L  —    Z-  X  o,  +  -2 ! 

a 

.   .  ,.    10 

tt  Jt  — ^  * 
4  t 

Vith  division         £>    ■-    -lection  io  fee***      Mmfcalnlai  #, 

and  substitution  of  0  m  n/V  ,    ifttPl  AT    is   c  ^naidere  itant, 

aation  10  reduces  to 
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3     +    S,« 

4        O  • 


J?  id-. 

d  .    11 


Prow  the   «ibov<  tnifclOBl   8   and  9 
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*■-  to       ' 

but  since  Jt  o  =  V  or 


je 
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In  th     a  asm  manner 
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Jtn  m  /,{!  «   L%2)   ■ 


o       ** v  .£    or 
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ITlth  these  and  division  bj/tT,    equation   11 

become  3 
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Prom  Equations   3,   0,   and  13  It  follows  that 
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3    — —    -   X  c 
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to  res  iunticn  lp  and  with  substitution  for 

A^  ^,    Equation  Ik  becomes 
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with  rearrangement 
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-  »  X    T    d  d 

>«>      0  Z  __        i  ■imi  iiiii  l      .*■  ■■■■»— 

'  *-i      rL  *■  d  fc 

,rco  f or eing,  k  will  be  a  eo.  uae  the 

pile  geometry  is  not   e  >4.  B   cori 

value  of  noutron  density  will  be 

n  ■  r^  ♦  n1e  J  to  i#   19 

where  n^  is  the  steady  component  and  n-.    is   the  complex  ampli- 
tude,     tm    in  analogc  nner 
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By  combination  of  ^.uations   .,  I  ,  and  2  ,   uation  lo 


4  /i  >.*VJ"*J 


(1  -  k) 


J- 


ly  ont   average  of  delayed  neutrons  will    >o   con- 

sidered in  -    to  ma -re   the   equation  more  wield} .      fhla 

introduces   so  clonal,   but  ■ via  11,    error  into   the  results 

(6,   p.   25D-      Loui^ion  21   in   the  complete  form,    out  by  the 
definition  of   a  transfer  function  only   t>ie   sinusoidal  portion 
is   of  conoorn.      -  1th  separation  of  t  dy  and  sinusoidal 

co-Tponents  and   with  the  station  dropped,   the 

sinusoidal  portion  of  Equation  cones 

jut  i-  Q_  j«*»t  -      T  -*  jot 


ju»  t  .22 

I 

(1   -    :■:)      ' 
atitutlon  of   Equations   19   and  20  into  Equation  15   give* 


1- 


*»i 


11 

The   sinusoidal  portion  of  liquation  ?3  beoones 


or 


fill  nx 


^fie-bT 


etltution  of   fixation  2$  irv  ;>ivea 

U>  t  V  J  «*>  t 


!iiO 


(k  -   1)  -   1) 


Lt>*j-»o,J 


\%  -  1)  (1  -  k) 

With  the  trr»  function  defined   *fl  tno  ratio  of  the 

response  to  the   ginusoidally  varying  driving  forco,    the 
desired  ratio   is 

nl 


1  -  k 
By  rearranf.  boa  26  :®s 

Tho  denominator  of  the   right)  side  oi"  -.tion  27  raay  be 

■iv  Led  to   girt   the  desired  final  form, 
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Figure   1  is   an  overall  view  of  the    «tubcritical  assembly 
the   east   cover  renovod.        o act or   grata    graphite 
cylinder 8  60  inches  and   7   inches   in  diaretor  were 

fir  chined  to  a  square  Torn  6   inches   across,    leaving  rounded 
oomore.     These  were   stacked   in  ten  columns  nine  rows  high. 
Additional  cyliriders  were  out  to   a  £-inc  by  6- inch  wide 

form  for  thf>  top   five  rows,      i'he  final    >ile   eift+j    with  the 
machined  cylinders   stacked  fourteen  rows  high  and  ten  columns 
wide,  treasures  60  by  60  by  79  Inches.     The  whole  array  of 
graphite   is    surrounded  by  a  »m  h  of   plywood,   0.  inch 

cad-  iu<-<,   ar.d  masonite,  bolted  in  sections   and  braced  with  a 
framework.      The  east  and  west  covers  are  removable   to  gain 
access  to  the   pile. 

The  Its   la   placed  on   a  wooden  oaee  constructed  so 

that  there   are  three  aoeessiole  spaces  ma  ;he  pile*      *.'nea« 

spaces   extend   the  north-south  width  of   the   assembly.     f| 
■1ml rtf  tu        ,    also  extending  the   width  of   the   pile,   were 
filled  with  water   ard  placed  in  the   two  outermost  apaoea  to 
moderate  the  neutrons    and  to  T4  radiations   coming  di- 

rectly from  the   sources  beneath  the   pile.  antral  space 

was  used  to  hou^e  the  oscillator  shown    In  figure  2. 
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Figure  1.  ?ha  suboritlcal  assembly 
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insure  2.      Thm  oscillator   assembly 
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Th»  boles  between  the  graphite  bloc      ,      '•asuring 
ately  1     Inchea   across,   pomlt   Insertion  of   the 
»ntfl   in  various    lattice   arrangement-* .      fclwilWf 
cylinders  appr  xi ■•lacel.y   C2   Lnehea   Ion;-;  each  containing  seven 
0,2-inc-i   uranium   slugs   constitute   the  fuel  elements* 

The   east  cover  was  in  ;  the   test  runs,   and 

five  ho  lea  were   (trilled  on  the  cover  center  lino  -tnoes 

of  one   foot,    tvo  feet,    three  feet,   four  feet,    and  five  feet 
from  the    oase.      For   all  further  discuseion  these  holes  will 
be  designated  by  ordinals   with   the  hole   one   foot  above  the 
base  call.  •   first   hole,    the   hoi*     two   feet   above 

the  base   being  called  thi  hole,   arid  so  on,  &h  the 

8*£~i  'attics,   as  shown  in  -  •   1,   th©  cover   holes  were 

in   line    with  empty  b&fce  bloclrs,   per- 

il neutron  detector  to  be   inserted  into  the  pile* 

The   apparr,  ;-<ed  was   assembled  from  readily  available 

material   and  parts.      fertrme  accuracy   and  refinement  of 
deni-.  n  were   not   atte?-    -      .  ve    Ptt-  e  soirees  cylindrical 

in  shape,   each  measuring  1  inch  in  &  '^er   and   1  3/8  inches 

high,   uQi-o   p laced   in  the   source  retainer   of   the  oscillator 
shown  in  figures  2  and  3.  -tlon  was   provided  by  a 

paraffin   jacket  formed  about  the  souree  retainer*      The 
jacket  was  made  by  rolling   l/l6~inch   alu  eet  to   an 

outside  diameter   of   7  3/8   inches  in;;   it   to   1/3- inch 

end  plates.         da  container  was  filled  by  means   of  a  filler 


V 

hole  with  25  a  of  hifrh  molting  point  v.  paraffin  and 

slowly  cooled  to  minimize  the  possibility  of  shrinkage  voids. 

rcximately  three   inches   of   paraffin  ounded  the   sources. 

The  retain© r  acted  as  suppoi  t  for    both  the   non-rotating 
paraffin  Jaciret   and  the  bearings  on  whic  external 

cylinder  rotate    . 

The  rotating  cylinder  was  formed  by  rol'  -inch 

aluminum  sheet   to  an  cut  aide  dimeter  of   8   ttM>h*S   and  cutt. 
It  to  a  length  of  1$  s.  ites,    also  of 

aluminum,   were  from    |    inch  stoc.,      hrln     fitted  crew< 

to  a  collar  fitted  on  the  bearing  outer  run  casing.     The 
cylinder  was   screwed  onto   Ifet   end   plates 

assembly.  ■  sine   pattern  was  cut   from  a  O.GlO-ineh  thick 

sheet  of  cadmium  to  the   dlmennions    shown  in  re  3   arid  oe- 

wented  with   an  aluminum  based  cement  to  the    inside  of  the 
rotating  cylinder. 

-'urther  discussion  position  1  will  be  used  to  refer 
to  the  rotating  ey  r   and   pattern  position  whenever  point 

A  of   the  flattened  pattern  in  i?i{?ure   3  is  top  dead  center. 
Position   2  will    jo  used  to   1  to  that  point  B   Is  t,  id 

center,      t'ased  on  the   difference   of  readings  with  the 

am  pattern  in  the   two  extreme  positions.    It   is   estimated 
t   approximately  30  it   of  the    source  neutrons   reached 

thermal  speeds   in  the  paraffin.      It   is   also  estimated  from 
the   accepted  macroscopic   absorption  cross   section  of  cadmium 
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Figure  3»   Oscillator  and  cadmium  pattern  detail 
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that  98  per  cent  of  the  thermal  neutrons  reaching  the  pattern 
were  absorbed* 

An  extension  of  one  of  the  bearinr  collars  was  -Trooved 
to  hold  the  leather  drivir     1 t.   The  rotation  frequency  was 
varied  by  controlling  the  volta.re  to  a  ;-lnch  Black  and  Decker 
electric  drill  mounted  on  a  four-foot  wooden  extension.  A 
speed  reducer  was  Interposed  to  permit  eneeds  down  to  approxi- 
mately 5  Rf  .   The  oscillator  was  desired  for  speeds  as  high 
as  1009  R?M#  but  no  testa  were  carried  out  in  this  range. 

Por  maximum  signal  the  neutron  detector,  a  3/3 -inch 
diameter,  8-inch  lon/$  tube  filled  with  enriched  ftp-  gas,  was 
inserted  in  the  first  hole  to  tie  oile  midpoint.  A  Radiation 
Instrument  Development  Laboratories,  model  200-1,  scaler  was 
used  for  the  hi#h  voltage  and  initial  amplification.  The 
amplifier-discriminator  jumper  was  disconnected  and  the 
amplifier  output  taken  directly  to  a  Nuclear  Chicago  model 
1615  survey  meter  with  the  high  voltage  rendered  lnooerative. 
The  output  of  the  survey  meter  was  fed  to  a  type  112  direct 
coupled  Tektronix  amplifier  and  thence  to  one  amplifier  and 
pen   of  a  two  pen  T3rush  recording  system.  The  instrument 
arrangement  is  pictured  In  figure  if,  and  all  the  circuitry 
utilized  is  shown  by  block  diagram  in  Figure  $. 

The  time  applied  variation  in  source  strength  was  re- 
corded by  means  of  brush  contacts  which  were  mounted  on  the 
oscillator  frame  and  the  rotor  and  which  were  wired  parallel 
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to   a  10 I  real Bt or   as   shown  In  resistor 

vaa   In  aerien   with  a    1001     res!    X      ,    and   a   110-volt   60- 
cycle  voltage  was  >ssed  across   ta<     .  o   seoond   ra- 

ce amplifier  ind  pen  were   also  connected   acrons  the    10) 

resistor,  .he   rotor   i  oto   ojon,    a   ID-volt  60-cycle 

>ut    signal  v  as   plotted  by   I  esu      Whenever   the   rotor 

co-  closed,    this  was  shorted  out.         m   oruah 

contact  nounted   on  the  fraene   was   approximately   120  degrees 

ot.  rotor   top   dead  center  vhiie   the  rotating  contact  was 
set    riMraftl   li   point   A  on  the  rotor.      It  was  found 
oscillator   speeds  up  to   about  300  could  be  treasured  in 

this  manner.      *t  faster   speeds   the   1  .i   of   tine    tar   wave 

is  shorted  is   only  a  portion  of   one  wave  cycle   and   is   very 
difficult   to  read   vita   t.o   rv;,her  slow  pen  responne. 

idings  were   taken  with   the   oscillator  steady   in  posi- 
tions   1    and  2  to   determine   the  range,      tttli    ninulated 
an  infinitely  alow  oscillation.      The   difference   in  the 
readings  was   called   S-^'.      a«adi--          ore  th€?n  taken  with   I 
oscillator  rotating  at  various    speeds,    and  for  each  fre- 
quency the  -peak  to  pea];  amplitude,    called  n^1,   was  compared 
with   the    oaso  value    •%*« 

The   phase   ancle  was    determined  by  plotti:  alt  Ion   1 

of  the  rotor   on  the   time  wave  and  comparing   it   with  t 
position  of  peat:  flux  for  each  cyole.      It  was  noticed  during 
preliminary  tents   that  moving  ttM    p*  ;eadlly  in  and  out 
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e   pilft   by  ha  ced  a  recorded  wave  having  an 

amplitude   less  taan  -iy  H  >© 

were  bfctfl  to   stop   In  Mat   pile  foj  to    ion  3ecor,d3. 

The   i>en,   hovev.  ,    ■#*   »4  to  reaoond  in  ;e  with  bho   action. 

was  therefore  dec;.. led  to   correct  she   recorded   anp  Liiu 
for   llM   two   second    tirae   c  survey   Kieber  and   to 

aooe  >t   Mm  recorded  phase  angle   as   oeing  approximately 
correct. 

'cii  circuit  analysis   &p>  eared  to 
too    Lengthy,    it.   was   dec  id,  correction 

en tally.      A  tint  wave   signal  rat  or  was   congee tad 

to   tilt?*   •  '  net*?      '.  id  jus  ted    to   give   a  ^-volt   out- 

put,     ^tie  scaler  ainpiir  as  not   included  in  the   teat   cir- 

cuit LSS  of  the   internal  circuit  >difieatlcr  it 

would  have   been   entailed  with  t  Lusion.       .dace  tne 

scalor  awpli/ ■■  in  essence  ti  voltages 

with   a  five  nicrosecond  reaolvbion,    it  was  iered  to  have 

negligible   effect  on  the  variation  of  the    i  ooorded  signal 
amplitude.         Mi    signal  gilUftHi   frequency  was   increased 
until  the   survey  meter  rag  od  20,000   counts  tnute 

which  was  approximate  i®  same  as  the  reading  outainea  with 

the   probe   inserted  to  the    rill  int    in  Mm   first  hole* 

It   Had   been  noted  as   the   oscillator  waa   P0t*t*4 

very  slowly,   the  survey  meter  rea    '        varie  approxi- 

mately  19,000  eour.  r  minute   and   21,00      i    ants   per  minute. 


The  signal  generator  utxu  adj,.  until  trio   survey  Hit 

read  values,    and   the   extremes  were  d  on  l&t   dial. 

The   sinusoidal  variation   jetwee;.  taes©  extr-  od 

by  ig  a  yardstick  aol.:  or 

dial   to   permit   a   pendulum   action,       i  li©avy  vol  olar  ped 

at  v;*  ,anees   on  the  >zr&  and   11  io  si;  va  late 

the   different   ayellc  ;ies.  rieid   in 

the   extreme  positions    Mm   Mplito  dfl   dlffereneej   as  aefore, 
air                 m   infinitely  slow  cycling   speod.      The   I  I  .io  of 

served  pea  any  o  to  the 

■■plittu  rence  d-  ibove  war-    I  I        MP* 

roction  faster   for-  that   :  ■■    .      Tho  a'~>ilr;ude   correction 

curve   la   plotted   in       .        o  6.  a  for   the 
on  Oirvo   imi    :r?v>;r-eXy   lir"  th. 

.  jgytns^pMB  ,  less 

inches  were  not  fi  Le*  net    I  us 

handle  [  ided  to 

.iter   cyoli:-  <eds  vhich  were  beyond  the  u.  ues  of 

I   •  a.      Hm   30    inch  limitation  of  the   yard- 

>     -     •  '  Mi   dial*   cc rrss  -ended   to   a 

radians  per  second  usir  .pie  pendul 

>-ory.      It   is  I  v,ed  that   to  obtain  a  ft?  of  one 

radia*  &  32»2-foot  pel  ri 

M   correction   mm  fti  by  dividing  nx»   by  the 
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Figure  6,     Amplitude  ratio  correction  factor  vs.   frequency 
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corresponding  correction  faotor    at   the  same  frequency  to  give 

arlmentul  values  of   amplitude  ratio,   nVo,1,    plotted 
in  Figure  7. 

Sotae   to  \s  were  utter  the  second  and  third 

hole*  but   the   data  wore  not  readable,      apparently  the 
statistical  variations  were  of  the  same  order  of  ma^nltu 
as  the  flux  variations,   masking   Hit   desired  signal* 
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The  primary  results  of  the  experiment   are  presented 
graphically  ll  iU'es  7   and  8.     'HhQ  data  from  vhic  se 

figure  a  were  drawn  are  tabulated  in  the   Appendix. 

*Rje   smooth  curves    ire  plots   of   Equation  25  for  different 
values  of   the  wan  neutron  lifetime   In  a  finite   assembly  and 
for   different  values  of  k.      olnce   neither  X,  nor  k  have  been 

;ermined  for  the   assembly  In   question,   values  of  .7 

and  G.B  and   X»  ■  0.07  seconds  and  0.1  seconds  were  chosen  to 
indicate  the   trer-.d  of  the  curves  with  both  JL  and  k  and  to 
indicate  th«  general  curve  shape.      In  trie  theoretical  deriva- 
tion the   six  groups   of  delayed  neutrons   were  replaced  with  an 
average  group  to  reduce  the   amount   of  computation,      v.lth  the 
total  fraction   p   of  the  delayed  neutrons;   equal  to  0.0073,    an 
aver  a  e    \   i<j**l  to  0.08   seconds"*     waa  computed  from 

\  -      P« — tt-  16,   p.   zyd) ,  th  the   aoove  eters 

in  the  liquation  28j   the    amplitude  ratio  and  phase  angle   were 
determined  for  various   frequencies. 

tfoasurements  of  tLe  neutron  flux  in  the  first  hole  were 
taken  for  oscillator   speeds  from  approximately  0«$   to  5 
radians   \mT  second.      ;  otn  phase   an^le   and  amplitude  ratio 
readings  were  obtained.      The    *.*ecorded  phase  angle  ia  plotted 
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Figure  7.     Amplitude  ratio  vs.   frequency 
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in   fj  pU*«  8.      The   recorded  amplitude  ratios   were   corrected 

tvo  see*  •    eCAftfcant   of  the   su;  ox-  by  the 

correction  factor  Indicated   in  &  6.      The   data  Tor   the 

•■plitude  ratio   correction  factor,    tabulated  in  tho  -ndix, 

plots  as  a  ■  tar  eight  line  on  "<ni- logarithmic  paper.       Itt 
range  was   increased  by  straight  line   extrapolation.  la 

result   is   valid   only  over  a  limited  •  since  for  all  fre- 

quencies the   factor  cannot   ee    grtStter   trim  unity  or   less   than 
zero,      hence  the  true  factor  must  deviate   somewhat  from  a 
straight    line   as  it    approaches  tho  extreme  values.      At  the 
lower   frequencies    the   correction  factor    itself   ;>ecomea  small, 
so  the   a^ove   deviation  will  have   little   effect  on   the   cor- 
rected amplitude  ratio.     At  tne  higher  frequencies  w:\ere  the 

rection  factor  becomes   large   In  comparison  to   the   readings 
the  effect   Is  magnified,   and   extreme  oa, e  would  need  to  be 
taken  to   justify  any  extrapolated  values,      fortunately  the 
fre  :Ies  obtained  in  the  exp  ntal  correction  factor 

tests  wore   as  high  as  those  used  in  the   amplitude  ratio  tests 
and  no  extrapolated  values   in  the  frequency   range  nad  to 

be  .      The  u)    plotted   in  Figure  6   is   actually  that  of  the 

jdulum,   but   it   corresponds  directly  to   that  of  the 
oscillator-   since   in  both   cases   it   is    the  frequency   of  the 
drivinr  function.      The   frequency   was   solved  for  a  simple 


pendulum   (without   damping)   fross  tho  formula    U)     =      lg/ 


■X 


:> 


where  f   is   tae   gravitational  acceleration   in  feet    per   second^ 
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Figure  8»      Phase   angle  vs«   frequency 


and    X-    Is   tie  pendulum   I  t«    • 

&B  into renting  result  of  the  c  -sent  was  the  decrease 

in  nedian  flu.x  level  with   Increaal     ;    CO  .     The  value  of  the 
median   line   of  the   sinusoidal  fix;  varieties  for  each  to 

frequency  was   plotted  in  figure  9   as  it  appeared  on  the  pen 
recorder,   with  tne  reco  >&per  scale   as  the  arbitrary 

ordinate  seal    .         .Is  value   correfc  to  n     in  the  theo- 

retical development*    and    it   is    tabulated   in  ondix  for 

each   fro  'U-ncy  tested*      In  the   split  bin,..:  of   Equation  21  into 
the  steady  and  the  sinusoid-il   condition-,    It   is    seen  that  n{ 
has  no  theoretical  rolati  CO.      lease  I  nervation 

was  not   predicted  by  -ory.  ioal  explana- 

tion lies  in  th©   oscillator  d— ign   and   t  ihod  chosen  to 

ncc  "3b    sour  of  forcing,      lio   tests  were  run  to  de  ternlne 

the   exact   c^une  I   the   actual  B  flux   IovjI  would   nave 

little   effect  on  I  cling  amplitude  ratio. 

!  jerirnental  uncertainty  results   f ro! I  A 

very  slow  rotation   trial  run  Ml  it   with  the  neutron  probe 

vltr.ii?   three   inches  of  the  oscillator  to  check   th©  variation 
?.n  the    SGttFfte*      Use  vox-  I   rtions   noted   lr. 

»tted   sine  wave  could  possibly  have    oeen  due   to   neutron 
scattering,  voids   in  the   paraffin,   and  fc. 

pS    the  greatest   source  of   orro.     In  the  eip  ST  lies  nil  was 
the   signal  f  I  *t   to   si  H  variations   and 

noise.       :-iese  variations   obscured  distorted    ,  oaks   of 
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Figure  9.      Median  flux   level  vs.   frequency 
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•are  10.      Views  of  typical  data 
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A.      Low  gain 


B.   High  gain 


ale,   adding  I  detei         Lug  both  the 

amp-  .i»o  angle,  o  samples   of  t^pioal 

data,   df  ri  [        readings  61  I   low 

gain,    are  shown  li  tra   10. 

,al   and   9  .od  values   if 

;.Io,    It   Is   or  io^s:ii7   to 

ooaervo  the  ■  '   ',    >   carvoa.      It   Is   saan  *« 

•x  fcal  polntn   Ho   reasonably  close  to  tUe»    theorek'cal 

vo  plot  tod  fo  i  of  ■   -      .      md   A-    -  ';,0?  da« 


The   transfer   fu  n  dote; -mined   ey  ontally  agrees 

its  veil  with   t  ^ed  by  the   the  or.  .       hile   it  was 

not   the   purpose   of  the  ,    the  exp«     i     r;tal  points   indi- 

cate a  mean  neutron  lifetime  in  the  assen'uly  of 
0.07  seconds   and  an  effective  multiplication  factor  of 
approximately  0.8. 

Hie  rotor   assemoly  t  Used  and  absorbed  a  sufficient 

percentage  of  neutrons   to  permit   obtaining  acceptable   data 
ft/mm  the   first  bale*      Si  second  an<  &v  holfs  the 

sinusoidal  variation  ir  neutron  flux  was  reduced  to  ttee 

order  of  magnitude  of   fchs   rt  ii   and  noise  variatic 

and  could  not  be  read.      I  life   •   larger  or  more  efficient 
neutron  detector,    a  thicker  paraffin  j  ,   or  possibly  a 

thicker  Sftdmltti   p*fct«r&j    the    amplitude  of  tie   sinusoidal 
variation?   could  oe  substantially  increased  permitting 
correspondingly  smaller  deviations   and  statistical  errors. 
Han   level   of  >ile   neutron  flux  decreased 

h  increasing   cyeli  need.      This  result   is  not  cted 

by  the   theory  and   is   p,    ■   ■    ly   due  to  the  iktf   the 

oscillator. 

exp©  tal  result;:  t  not  extremely 


accurate   It    la  felt   that  the   data  can  be  Lth 

*uf  -'fccy  and   go  iy  to  tantiate   the 

theoretical  predictions. 


Ul 


I 

It   Is  very  apparent   that   the  possibilities   In   eye I I 
a  suberi  assembly  have  not  boon  fully  exploited.      It  re- 

mains  to  improve  on  and  enlarge  t  of  reading   to  nuoh 

fas  to,  much  .  is  will  require 

some  modification  In  the  oscillator  drive  -lechania :  . 

.   a  higher  degree  of  certain  the   greater 

speeds   can  bo   c  ning  iftftifl   circuitry  to 

have  a  lower   time   c  to  Include   a  di fi- 

eri or  to  r»4u««  iml  rarl  .         -ce  this 

has  been  accompli'  It  v;cald   be  desirable  t  sJLa  the 

sp*  llo«  of  tiw  |  i*  of  re- 

activit  ting   and  to  compare  the   results  of  tae  tv 

r.ethods* 

HI   t-K:  --action.*-    aaove  certain  to  the  assembly 

itself   sj    !!oing  the  object  of   Investigation.      MS  e 

background   knowledge  of   the    transfer  function,    the  whole 
area  of   tnvosl  :  the  utilisation  of  the   subcriticil 

assembly  as  a  facility  fo  iting  c  <*nts   and  for 

advanced  c  rese  will   have  been  opened. 


i.      illardlta,  Ccrbi    ,  til, 

a.  Commlaaion     [Ttoou— nt] 

2      [  ]  •  ;®» 

on, i      .       »  EH  I         Com-nlasion.        M«    13» 

2.  ixtnaan,      .      .,    Ifeaaauer,    G. , 

activity  meaaurom*ntfi   In  a  subcritloal  pile. 

ommiaaion        [  l]  -  [  . 

I.   6n  ■  H amour 9  and        .]    .  uata,   Ga., 

iver   I^aboratc  ■  7.  j.    1S>55» 

3.  -el:.     ',       ...  .      .,  ,  casura- 

nfc  of   the   transfer  function  of  a  fast   critical 

Laboratory      [  nt] 

2.  .or,       .       .,    uITioe   of    f<  cal 

vices. 

.      ^rovn,    border    •,.  »11«  .  Lfti  of 

amrvone  .  S*  X*i    John     iley  and 

,    I  ?:c .  I r8 . 

£•  ,    Joseph    P.  3.      a  tor  ic 

>sion  tea  if  lad      [Document] 

3260.      Oak     tldgtj  -•■.'-. ,    .'eohnloal   Information 

;rvice#      July  16,  1  .,  v. 

6,      Glass tone,    ta uel.      ?rlncir>lea  of  nuclear  reactor 

englneerir}.  .        H    .'ton,    <:.    J.,     ..    .         km  triad 
Co.,    Inc.      1955  « 

7-     -_-_»____  ilund,    Hilton      .  M&f   of  nuclear 

ictor  theory  .        .  ineeton,   M«    J.,      .  -.'.rand 

Co.,    Inc.      1 

8.      Barrar«    J.    .'.. ,   hoy**',      .      . ,    and   r.rucoff,    D.  Ter 

function  oi  -2  roactor.  i loonies    iOj 

no.   8s   32. 


1*3 


9.      Oimag    J.    I.,    Jr.,      .  "vc    ,      .      . .  and   i  igott,    •-•. 

©d  auto  r»tera  for   b 

ectady  raactv     .       itonilo  .-nergy  Commission  D«- 

isslflfri     [                   ]  -                            idg», 

.,  ^tension.      Haroh  29 i 


J-V 


9. 


10.  . ,  Mm  ,     .   I.     «  it 

or  plant*  --.iralsslo::     assif iod 

[        ]       -11,'         , 
Technical  Information  .service*     y  11,  19§1. 

11.  ,   .   .   Control  of  nuclear  reactors  and  power 

>5. 


hk 


The      r  wishes  to  thank  Dr.  Robr  I   .  TJhrl^  for  the 
ordinal  ■mgytl—  of  tfcn     '.].?ct  und   for  his  interest  and 
patient  assistance,   any  thanks  are  due  Or.  Glenn  Murphy, 
Head  of  the  depart  lent  of  Thsoyti— 1  &nd   Applied  Heohanios 
and  Chairman  of  the  Uuclear  !ngineerlM  administrative 
Cownittee,  for  his  assistance  in  securing  the  instrumentation 
required  for  1    etion  of  this  work.   Aa  is  true  in  any 
project  undertaken,  the  associations  and  occasional  meetings 
are  too  numerous  for  an  individual  recountal,  but  an  acknowl- 
edgment would  not  be  complete  without  expressing  sincere 
gratitude  to  t  ie  many  persons  who  have  given  suggestions  and 
encouragement  throughout  the  entire  course  of  the  work. 

The  author  would  also  like  to  express  his  appreciation 
to  the  0*  S.  Dfaval  Postgraduate  'School  for  the  educational 
opportunities  afforded  and   to  the  United  states  Mavy  for 
making  all  fchio  work  possible. 


us 


Table   1.     Amplitude  ratio  action   factors 


1  cS  U>  rreotl 

Inches     feet     (=  %/l)     radians/second     amplitu 


cfcion 
factor 


30 

.5 

12*57 

3.59 

.95 

0., 

25 

,08 

15.5 

3.9k 

0.75 

20 

!••' 

19.35 

k* 

0.6 

.12 

15 

1. 

.7 

5.08 

0.35 

.07 

Table  2.      Aasplltud©  ratio  data 


him 


ians/r»ocond 


V 


n, 


,  • 


4tf 

0. 

0.85 

.95 

13 

0.87 

.65 

0.65 

.92 

Ik 

1.- 

0.5 

0.5 

0.88 

15 

1.73 

O.i* 

0.1* 

0. 

16 

2.25 

0.3 

0.3 

0.1 

17 

.78 

G. 

0. 

.69 

18 

3.31 

.15 

0.15 

0.65 

19 

3.73 

0.1 

0.1 

.59 

Tabla  3.   Phase  angle  data 


Hun  tj  Phase  an r.le -radians 

radi ana/ 8 eoond 


o.m*.  .07 

0.73  0.27 

1.15  O.kk 

U7J  0.53 

1.65  0.62 

3.31  0.38 


fable  1+.   Median  flux  levels 


co  Level 

radiems/aecond 


0.1*9  7.3 

0,37  7.3 

1.26  7.25 

1.73  7.2 

2.70  7.2 

3.31  7.15 

3.73  7.3 

%  6»9 

98  7.0 

9.67  7.0 

12.2  7.0 

13.9  6.9 

16.1*  6.8 

17.9  6.3 


/ 


tl 
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